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We report on a search for the flavor-changing neutral current decay fJ'^f^^ 
in pp collisions at ^/s=1.96 TeV using 65pb~^ of data collected by the CDF II 
experiment at the Fermilab Tevatron Collider. A displaced-track trigger se- 
lects long-lived candidates in the — » search channel, the kine- 
matically similar tt^tt^ channel used for normalization, the Cabbibo- 
favored —>■ K^n^ channel used to optimize the selection criteria in an un- 
biased manner, and their charge conjugates. Finding no signal events in the 
—>■ fi'^fJ^~ search window, we set an upper limit on the branching fraction 
B{D^ fi+fi-) < 2.5 X 10"*^ (3.3 x IQ-^) at the 90% (95%) confidence level. 

PACS numbers: 13.20.Fc, 14.40.Lb 

The flavor-changing neutral current (FCNC) decay D'^ f^^lJ'~ ti] is highly suppressed in 
the Standard Model (SM) by the nearly exact Glashow-Iliopoulos-Maiani (GIM) 3| cancel- 
lation. Observation of this decay at a rate significantly exceeding the SM expectation would 
indicate the presence of non-SM particles or couplings. In the context of the SM, Burdman 
et al. B] calculate the branching fraction to be B{D^ —>■ fj^^fi^) ~ 10^^^ from short- distance 
processes, increasing to i3(i5° /W^/U ) ~ 10^^^ when long-distance processes are included. 
This prediction is many orders of magnitude beyond the reach of the present generation of 
experiments, whose most stringent published limits are 4.1 x 10~® from BEATRICE j^j and 
4.2 X 10- from E771 1^ at the 90% confidence level. Thus, a large, unexplored region exists 
in which to search for new physics. 

Burdman et al. consider the effects on fi'^fi' from a number of extensions to the 

Standard Model: R-parity violating SUSY, multiple Higgs doublets, extra fermions, extra 
dimensions, and extended technicolor. They find that the — > fi~^fj^~ branching ratio can 
be enhanced by orders of magnitude to the range of 10~^ to 10"^'' in these scenarios, and in 
the case of R-parity violating SUSY, roughly to the level of the existing experimental limit. 
Similar enhancements can occur in K and -B-decays, but charm decays provide a unique 
laboratory to search for new physics couplings in the up-quark sector. 



This search uses a 65pb~^ data sample recorded by the upgraded Colhder Detector at 
Fermilab (CDF II) at the Tevatron pp collider with y/s = 1.96 TeV between February 2002 
and January 2003. The components of the CDF II detector pertinent to this analysis are 
described briefly below. Detailed descriptions can be found elsewhere CDF uses a 
cylindrical coordinate system in which is the azimuthal angle, r is the radius from the 
nominal beamline, and z points in the proton beam direction and is zero at the center of the 
detector. The transverse plane is the plane perpendicular to the z axis. The pseudorapidity 
ri is deflned as r? = tanh~^(cos where 9 is the polar angle measured from the z axis. A 
silicon microstrip detector (SVX II) iTf] and a cylindrical drift chamber (COT) [S| immersed 
in a 1.4 T solenoidal magnetic fleld track charged particles in the range |?7| < 1.0. The 
SVX II provides up to flve r-0 position measurements, each of roughly 15 fim precision, 
at radii between 2.5 and 10.6 cm. The COT has 96 measurement layers, between 40 cm 
and 137 cm in radius, organized into alternating axial and ±2° stereo superlayers. The 
solenoid covers r < 150 cm, and electromagnetic and hadronic calorimetry occupy the 
region between 150 and 350 cm in radius. Four layers of planar drift chambers (CMU) ^ 
outside the hadron calorimeter cover the region \ri\ < 0.6 and detect muons of transverse 
momentum > 1.4GeV/c penetrating the 5 absorption lengths of calorimeter material. 

The decays used in this analysis are selected with a three-level trigger system. At 
the flrst level, charged tracks are reconstructed in the COT transverse plane by a hardware 
processor (XFT) |lO]. The trigger requires two oppositely charged tracks with reconstructed 
transverse momenta px > 2GeV/c and pti + Pt2 > 5.5GeV/c. At the second level, the 
Silicon Vertex Tracker (SVT) jll| associates SVX II position measurements with XFT tracks. 
The impact parameter of the track, d^, with respect to the beamline, is measured with 50 fim 
resolution, which includes a ~ 30 fim contribution from the transverse beam size. Requiring 
two tracks with 120 /im < |(io| < 1.0mm selects a sample enriched in heavy flavor. The two 
trigger tracks must have an opening angle satisfying 2° < | A0| < 90° and be consistent with 
the decay of a particle traveling a transverse distance L^y > 200 /im from the beamline. At 
the third level, a computing farm performs complete event reconstruction. The sample of 
~ 10^ D*-tagged two-body D° decays selected by the trigger is used to estimate backgrounds, 
to optimize selection requirements, and to normalize the sensitivity of the search from the 
data sample itself. 
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The D'^ — > branching ratio, or upper hmit, is determined using 

^ r A- y - V ; iV(7r7r) e(/i/i) a(/i/i) ^ ' 

where B(D^ tt+tt^) = (1.43 ± 0.07) x 10^'^ is the measured normahzation branching 
fraction jl^, N(fifi) and A^(7r7r) are the numbers of D° ^-iid D'^ tt+tt" events 

observed, and e and a are the efficiency and acceptance for each mode. Except for the 
requirement of muon identification, and the assignment of different particle masses, the 
same selection requirements are apphed to both modes. In this analysis, we determined 
the upper limit on the number of signal events observed, N{fifi), by assuming that the 
number of events found in the signal region is the sum of signal and background events, both 
obeying Poisson statistics. Normalization was made to D° — > vt+tt^ rather than the more 
numerous K~tt~^ decays. Kinematically, the D'^ tt+tt^ mode is nearly identical 

to A*^/^") minimizing the differences in acceptance and efficiency, and introducing 

minimal systematic uncertainty to the result. The width of the reconstructed mass peak for 
two-body decays of the in CDF 11 is about 10 MeV/c^, sufficient to separate K^7r~^ 
kinematically from D'^ tt+tt^ (Fig. Q). 

In the spirit of obtaining an unbiased result, a "blinded" analysis was performed. The data 
in the signal mass window were hidden and the analysis cuts optimized without knowledge of 
their actual impact on the result. The optimization was performed on kinematically similar 
but statistically independent events. Only after all selection criteria had been fixed was the 
signal region "unblinded" and the final result determined. 

We first outline the general event selection requirements common to all the data samples 
used in the analysis and then discuss how they are used to determine the quantities in Eq. (pj. 
All of the samples consist of D*^ — > D^ti^ candidate decays coming from data sets where all 
requisite detector components were functioning properly, specifically, the SVX, COT, and 
CMU detectors, and the displaced-track trigger chain. candidates were formed from pairs 
of oppositely charged "trigger tracks" that are "CMU fiducial". Trigger tracks are tracks 
reconstructed offline that have been matched to online SVT tracks. A track that intercepts 
the active region of the CMU when extrapolated from the COT through the magnetic field 
of the detector is said to be CMU fiducial. 

To select D'^ candidates in a given decay mode, K~^7[~, 7i~^tt~, or we evaluated 

the invariant mass of each pair of trigger tracks using the corresponding mass assignment 
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and kept candidates in the range 1.840 GeV/c^ < Mpair < 1.884 GeV/c^. This corresponds 
to shghtly more than ±2cr around the mean of the mass peak, 1.862 GeV/c^. The D* 
tag reduces non-D^ backgrounds and ehminates the mass peak in the Kir channel due to 
mis- assignment of the K and vr masses. D*^ ^ D^^^t decays were selected by combining an 
additional pion track (tTs) with the candidate and requiring the mass difference Mpair+Tr^ ~ 
Mpair to lie in the range 144MeV/c^ to 147MeV/c^. The vr^ track was not required to be 
CMU fiducial or to be a trigger track, but it had to have the Cabbibo-favored charge for the 
Kti decay. 

The ratio e{7CTT)/e{^^) was determined from the muon identification efficiency and the 
pion reconstruction efficiency, measured in other analyses, as follows. From a sample of 
J/ip ^ fi~^fi~ decays collected by a trigger requiring one identified muon and one SVT 
track, the CMU identification efficiency for the unbiased muon was measured offline as a 
function of its transverse momentum px- We convoluted the efficiency spectrum with the 
Pt spectrum of pions from tt+tt^ and determined the effective dimuon identification 

efficiency to be e(yU/i) = 0.800±0.030. Using a detailed GEANT jl^ detector simulation, the 
pion reconstruction efficiency was found to be 95±1%, yielding e(7r7r) = 0.90±0.02, where the 
inefficiency arises primarily from hadronic interactions with detector material. Combining 
these values we find e(7r7r)/e(/i/i) = 1.13 ±0.04. Using the same detector simulation, we find 
the acceptance ratio a^nn) / a{fifi) = 0.96 ± 0.02. 

The number of tt+tt^ decays, N{7ctt), was determined by fitting the peak in the 

TTTT invariant mass spectrum. We performed a binned fif with Gaussian signal plus linear 
background, as shown in Fig. ^ Both the mean and width of the Gaussian were free 
parameters in the fit. N{inT) is the integral of the Gaussian over the ±22MeV/c^ mass 
window around 1.862 GeV/c^. 

The background to fi~^fi~ was taken as the sum of two contributions having different 

mass spectra: a peaked contribution from D° tt+tt^ decays in which both pions are 
misidentified as muons, and a relatively flat background due to all other sources. The 
flat background was estimated from the number of fifi candidates in a high mass sideband 
spanning the range 1.90 GeV/c^ < M^^ < 2.05 GeV/c^ with both tracks identified as muons. 
Before muon identification of the tracks is required, the distribution of events in the high 
mass sideband is found to be roughly constant, and we assume that this remains true after 
requiring muon identification. The expected fiat background is the number of sideband 
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events scaled by the ratio of the width of the signal region to the sideband region, 44/150. 

The misidentification background was estimated from the number of tt+tt^ events 

reconstructed with the fifi mass assignment and lying in the ±22MeV/c^ signal window 
(shaded area falling between the arrows in Fig. Q]) times the square of the probability for a 
pion to be misidentified as a muon. The vr-misidentification probability was determined from 
the sample of D*-tagged K~tt^ events. The average vr-misidentification probability 

is 1.3 ±0.1%. 

Three additional selection requirements were imposed. To remove instances in which the 
two D*^ decay daughters extrapolate to the same region of the CMU, potentially correlating 
the muon identification of the two tracks, we cut on the azimuthal angle AipcMU between 
their projections into the CMU. To suppress combinatoric backgrounds, we cut on the impact 
parameter with respect to the beamline, dxy, of the reconstructed trajectory. Further, 
we cut on the transverse decay length of the candidate, L^y The values of these cuts 
were optimized as described below. 

We determined the optimal cut values by maximizing a figure of merit given by 
S'/(1.5 + \fB) IJ] where S and B represent the number of signal and background events, 
respectively. This quantity has desirable properties for an analysis where both the signal 
and background are small: it behaves as S/ \/B for large B and it behaves as S as the esti- 
mated background approaches zero. The constant in the denominator is chosen to favor cuts 
that maximize the discovery reach at 3cr significance. To estimate S in the optimization, 
we used the tt+tt^ sample. To estimate the misidentification component of -B, we 

used a sample of K^-k^ decays in which both tracks were found to be misidentified as 

muons. To estimate all remaining contributions to B, we used the subset of the high-mass 
TTTT sideband sample in which one track was identified as a muon and the other was not. 
Note that the events used to estimate B in the optimization are distinct from the events 
used in the final background estimate for the result. The resulting selection requirements 
are: \/S.(t)cMu\ > 0.085 rad, \dxy\ < 150 /im, and L^y < 0.45 cm. When applied to the samples 
used for optimization these cuts remove approximately 58% of the background events and 
12% of the signal events. 

Using the optimized selection requirements, 5.0 ± 2.2 events remain in the high mass 
sideband, yielding 1.6 ± 0.7 expected from the fiat component of the background. The 
number of tt+tt^ events falling in the signal window when reconstructed using the 
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muon mass, multiplied by the square of the 1.3% vr-misidentification probability, yields 
0.22 ± 0.02 expected misidentification events. The total expected background is 1.8 ± 0.7 
events. The number of events in the normalization mode is N{TnT) = 1412 ± 54 (Fig. U) . 
Using this background estimate and normalization, the 90% confidence level sensitivity [l5| 
is 4.4 X 10-^. 

We apply the optimized selection requirements to the signal region of the /i/i sample 
and find no events remaining, as displayed in Fig. [21 Conservatively taking the number of 
background events to equal zero, the 90% (95%) confidence level upper limit on the number 
of —* fi^fi^ events is 2.3 (3.0). Using Eq. (P) we find an upper limit on the branching 
fraction of B{D^ n+fi-) < 2.5 x 10"^ (3.3 x 10"^) at the 90% (95%) confidence level. 

The uncertainties on N{7nT), e(7r7r)/e(/i/i), a{'n-7T)/a{fifi), and B{D^ —>■ tt^tt^) are incor- 
porated into the limit using the prescription of Cousins and Highland However, all of 
the uncertainties are smaller than 5% and have a negligible effect on the limit. 

In summary, we have searched for the FCNC decay f^~^f^~, using the new displaced- 

track trigger of the CDF II experiment. This is the first result from CDF in the field of 
rare charm decays. To minimize bias in the event selection, a blinded search was performed. 
To minimize dependence on Monte Carlo simulation, most of the needed quantities were 
determined directly from the data. No events were observed and we set an upper limit on 
the branching ratio of 

B{D^ fi+fi-) < 2.5 X 10-^(3.3 x 10"^) (2) 
at the 90% (95%) confidence level. This result improves on the best limits published to date. 
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FIG. 1: The mass distribution of candidate — > vr^vr" events. The ^^/^^ branching 

ratio was normahzed to the kinematically similar mode vr+vr^. The arrows indicate the 

±22 MeV/c^ mass window used for the signal. The curve is a fit over the range 1.82 to 1.90 GeV/c^ 
with Gaussian signal plus linear background. The shaded Gaussian represents the effect of recon- 
structing the events with a ji^ fi" mass assignment. The large Ktt signal below 1.82GeV/c^ is 
kinematically separate from the region of interest. The distribution of events in the region above 
the mass is roughly flat. 
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FIG. 2: The mass distribution of candidate — > /U^^^ events. No events remain in the 
mass region satisfying the event requirements. The events in the high mass sideband were used to 
estimate the background from all sources other than misidentification of vr^vr^. 
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